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Motivation 

Hypersonic  propulsion  system 


Ignition  time  (~10ms) 


»1 


Flow  residence  time  (~1ms) 


FI 35  engine:  (F35,  2011) 


Ignition  instability 


Challenges: 

•  Ignition  time,  Ignition  energy 

•  Flame  stabilization 

•  Combustion  completion 


Plasma  assisted  combustion 


Understanding:  Good 


poor 


marginal 
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Temperature 


Change  of  ignition  and  extinction  diagram:  the  S-curve  transition 
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Residence  time 

Scramjet 
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Research  goals 


Understand  the  fundamental  enhancement 
mechanism  of  plasma-flame  chemistry 

Develop  new  experimental  tools  to  validate 
plasma  flame  kinetic  mechanism 

Develop  numerical  methods  to  achieve  efficient 
modeling  of  detailed  plasma  flame  chemistry 


^5)  Outline 

1.  Background 

2.  Experimental  investigations 

•  Effects  of  plasma  assisted  fuel  oxidation  on 
flame  extinction 

•  Effects  of  in  situ  plasma  discharge  on 
ignition  enhancement 

•  Molecular  beam  mass  spectrometry  study  of 
low  temperature  chemistry 

3.  Conclusion  and  future  work 
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Background  and  previous  study:  flame  extinction 
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1.  Silicon  Controlled  Rectifier,  2.  Silicon  carbide  heater,  3.  R-type 
thermocouple,  4.  Fuel  injection  spacer  5.  MGA  plasma  power  supply,  5. 
MGA  device,  6.  MGA  power  supply,  7.  Cathode,  8.  Anode,  9.  Magnets, 
10.  Gliding  arc  initiation  wire,  11.  MGA,  12.  Insulator,  13.  Nozzle  with 
N2  co-flow,  14.  K-type  thermocouple  &  FT-IR  probe,  15.  Diffusion 
flame,  16.  Water-cooled  nozzle  with  N2  co-flow. 
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Only  thermal  effect! 

Ombrello,  et  al,  AIAA  J,  2006 


Previous  work  -  Ignition  study 


CtL/air  counterflow  diffusion  flame 


H2/air  counterflow  diffusion  flame 
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NO,  catalytic  effect 

1.  non  in  situ  discharge 

2.  Short  life  times  of  radicals  and  excites  species 


Ombrello,  et  al,  IEEE  Plasma  Sci,  2008 


Previous  researches  -  03 
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Mixture  fraction  gradient  dY  F  IdR 
(~  1/  axial  distance) 


Flame  speed  extraction 
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Ombrello,  et  al,  CNF,  2010 


Enhancement  [%] 


Concentration  [ppm] 


Previous 


researches  -  02(a1Ag) 


[°2(alAg)]>  ppm 

AHl,  mm 

3137 

4.76 

4470 

6.82 

4627 

6.83 

5098 

7.31 

02  (a^)  +  H  =  OH+O  fast 
02  +  H  =  OH  +0  slow 

~  5000  ppm  02(a1Ag)  4  2-3  %  Lifted  Flame  Speed  Enhancement 
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Change  of  Flame  Liftoff  Height,  AHL  [cm] 

Ombrello,  et  al,  CNF,  2010 


•  SDO  (w/NO) 
A  SDO  (w/o  NO) 
B  03  (w/o  NO) 


Energy  Coupling  Into  Flow 

~  1  eV  to  produce  02(a1Ag) 


10 


Previous  researches 


Crossover  T:  900  K 


O  quenched  even  at  60  Torr: 

How  to  utilize  radicals  efficiently? 


Sun,  et  al,  PCI,  2010 


Atomic  oxygen  effect 
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Pulse  repetition  frequency  (KHz) 
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Research  focus  in  the  second  year 


Thrust  1 .  Kinetic  effects  of  non-equilibrium  plasma-assisted  fuel 
oxidation  on  diffusion  flame  extinction  limits 

Thrust  2.  Direct  ignition  and  the  S-curve  transition  by  in  situ  nano¬ 
second  pulsed  discharge 

Thrust  3.  Plasma  flame  chemistry  study  in  a  flow  reactor  with 
Molecular  Beam  sampling  Mass  Spectrum  (MBMS) 

Thrust  4.  Development  of  a  plasma  assisted  jet  stirred  reactor  with 
molecular  beam  sampling  and  a  high  pressure  ignition 
chamber 
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Thrust  1.  Kinetic  effects  of  non-equilibrium  plasma-assisted  fuel 
oxidation  on  diffusion  flame  extinction  limits 
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Experimental  setup 


CK/Ar/He/CEL 


15  mm  x  22  mm 


10  mm 


10  mm  away  from  exit 


20  &  28  mm  ID 


U0  yfPo 


\ 


Time  (ns)  14 

The  thermocouple  was  coated  with  MgO  and  covered  with  grounded  Nickel-Chrome  sheath  to  remove  EMI 


Laser  diagnostics  schematic 
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Numerical  model 


Ar/He/02/CH4 

(0.32/0.4/0.26/0.02) 


discharge 
electrode 


Counterflow  nozzle  exit 
/ 


22  mm  10  mm 


Species  concentrations 
from  simulation 


Kinetic  model:  OSU  air 

plasma  model  [1 ,2]  with  USC 
mech  II  in  addition  of 
Ar/He/CH4  related  reactions. 

Physical  model:  quasi-one 
-dimensional  flow  equation  + 
steady  two-term  expansion 
Boltzmann  equation  [1] 


Reactions  [1-8] 

Rate  Const 

(cm3s'1) 

e+  02"^  e+20 

/( E/N) 

e+  02"^  e  +  O  +  O(D) 

/( E/N) 

e  +  CH4  CH3  +  H  +  e 

/(E/N) 

e  +  Ar  Ar*  +  e 

/(E/N) 

e  +  Ar-^  Ar(+)  +  2e 

/(E/N) 

e  +  He  He*  +e 

/(E/N) 

e  +  He  He(+)  +  2e 

/(E/N) 

Ar*  +  CH4  -»  Ar  +CH2  +2H 

3.3xlO-10 

Ar*  +  CH4  Ar  +CH  +H,  +  H 

5.8x10-1° 

Reactions 

Rate  Const 

(cmV1) 

Ar(+)  +  CH4  Ar  +CH3  (+)  +  H 

6.5x10-1° 

Ar(+)  +  CH4  Ar  +CH2  (+)  +  H2 

1.4x10-1° 

Ar*  +  CH4  Ar  +CH3  +  H 

5.8x10-1° 

Ar*  +  CH4  Ar  +CH2  +H2 

5.8xlO-10 

He(+)  +02  ->  0(+)  +  O  +  He 

0.6xl0-nT°-5 

Ar*  +  02-^  Ar+20 

2xl0-i° 

He(+)  +02(a)  0(+)  +  O  +  He 

0.6xl0-nT°-5 

He+20  He*  +  02 

1x10  33 

He*  +  CH4  CH  +  H2  +  H+  He 

5.6xl0-i3 

Reference: 

[1].  A.  Bao,  Ph.D  thesis  (2008)  OSU  [2].  M.  Uddi  et  al,  PCI  32(2009)  929  [3].  I.N.  Kosarov  et  al,  C&F  156(2009)  221  [4].  A.  Hicks  et  al,  JPD,  38(2005)  3812  [5].  D|_.g>. 
Stafford  et  al,  JAP,  96(2004)  2451  [6].  M.  Tsuji  et  al,  JCP,  94(1991)  277  [7].  A.M.  Starik  et  al,  C&F,  157(2010)  313  [8].  I.N.  Kosarev  et  al,  C&F  154(2008)  569 
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Experimental  observations  of  discharges 

02(0.26)/Ar(0.32)/He(0.4)/CH4(0.02) 


Absorbance 


Discharge  repetition  effect  on  species  concentrations 

FTIR  spectrum  with  different  pulse  frequency 
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Discharge  repetition  effect  on  species  concentrations 
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Carbon  deficiency:  5% 

Relative  uncertainties: 

<1%  for  CH4,  CO,  C02 
5%  for  H20  and  H2 

The  uncertainty  of  CH20  measurement  is  80  ppm 
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Reaction  path  analysis-CH4&H2 
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Reaction  path  analysis-H&O 


CH3  +  0 


OH  +  H 


CO  +  OH 


HCO  +  FLO/O 


CH4  +Ar*/e/Ar(+) 


3.2% 


O  +  H, 


Ar*  +  O 


02(a1Ag)/02(b1i:)  +  H 


+  ch4,  ho2, 
ch2o,h2, 

HCO 


OH 

H 

o2 

/=  40  kHz 
P  =  60  Torr 
T  =  300  to  933  K 


Mechanism  was  not  validated  below  700  K 
Large  uncertainty  at  low  temperature 

the  reaction  rate  at  300  K  for 
Of  D)  +  H2  =  H  +  OH  (4.4xl010  /cm3s) 
is  much  larger  than 
O  +  H2  =  H  +  OH  (2.6xl03  /cm3s).  21 


Extinction  limit  measurement  &  calculation 
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Faster  fuel  oxidization,  larger  extinction  extension 
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0.34  0.35  0.36  0.37  0.38  0.39  0.40 
Fuel  mole  fraction 


Simulations  were  performed  with 
experimentally  measured  boundary  conditions. 

OH,  H  concentrations  were  estimated  from 
simulation  by  matching  O  concentrations. 
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- simulation  (To=  423  K) 
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with  plasma  (/*=  10  kHz) 

- simulation  (To=  613  K) 

Oxidization  or  reforming  ratio 

Case  1:  fuel  was  oxidized  to  C02  &  H20 
Case  2:  fuel  was  reformed  to  CO  &  H2 

Fuel  reforming  enhancement:  fast  H2 
chemistry 

Fuel  oxidization  enhancement:  extracting 
chemical  enthalpy  rapidly 


Extinction  limit  measurement  &  calculation 


CH4  oxidization  ratio  (or  f)  increased ,  extinction  limits  increased  significantly 


5.3%  enhancement  from  H2 


The  dominant  enhancement 
mechanism  is  plasma  introduced 
rapid  fuel  oxidization. 

Deviation  is  due  to  additional 
reaction  paths,  but  not 
signif  icant  (10%). 


Simulations  were  performed  with 
experimentally  measured  boundary  conditions. 

OH,  H  concentrations  were  estimated  from 
simulation  by  matching  O  concentrations. 


Thrust  2.  Direct  ignition  and  the  S-curve  transition  by  in  situ  nano¬ 
second  pulsed  discharge 
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Experimental  setup 


02/He 


P 

f 


72  Torr 
24  kHz 


Laser  beam 


25.4  mm 


Power  ~  17  W 
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ICCD  images 


(a)  ICCD  image,  He/02  (0. 6:0.4)  and  He/CH4  (0.75:0.25),  50  ns  gate 

(b)  ICCD  image,  He/02  (0. 6:0.4)  and  He/CH4  (0.86:0.14),  50  ns  gate 

(c)  direct  photo  of  (a),  50  ms  exposure  time 

(d)  direct  photo  of  (b),  50  ms  exposure  time 
P  =  72  Torr,  f=  24  kHz,  a  =  175  1/s 
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OH*  emission  intensity  (a.u.) 


Classical  S-curve 


hysteresis  between  ignition  and  extinction:  S  curve 
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1400  ^ 
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1300  3 

Rayleigh  Scattering^12! 
1200  a.  method  for  T 
lioo  |  measurement  at  532 
1000  s  nm  from  Nd’.YAG  laser 
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900  *3 

800  S 

700 

0.05  0.10  0.15  0.20  0.25  0.30  0.35  0.40 

Fuel  mole  fraction 

Relationship  between  OH*  emission  intensity,  local  maximum  temperature  and  fuel  mole 
fraction,  T0=650  K,  Tf=600  K  He/02  =  0.66:0.34  ,  P  =  72  Torr,  f-  24  kHz,  a  =  400  1/s 


[1]  R.B.  Miles,  W.R.  Lempert,  J.N.  Forkey,  Meas.  Sci.  Technol.  2001  [2]J.A.  Sutton,  J.F.  Driscoll,  Exp  Fluids  2006 
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OH*  emission  intensity  (a.u.) 


S  curve  transition 


Relationship  between  OH*  emission  intensity,  local  maximum  temperature  and 
fuel  mole  f  raction,  P  =  72  Torr,  f-  24  kHz,  a  =  400  1/s 

He/02  =  045:0.55  He/02  =  0.38:0.62 


Fuel  mole  fraction  Fuel  mole  fraction 

monotonic  ignition  and  extinction  curve 
(monotonic5  curve) 


ignition  and  extinction  points  were 
pushed  to  lower  fuel  concentrations 
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Maximum  temperature  (K) 


Numerical  modeling 


OPPDIF  +  electron  impact 

Kinetic  mechanism:  USC  mech  II  +  OSU  air  plasma  model[1] 


E:  electric  field,  N:  particle  density 
Rate  constants:  Boltzmann  equation  solver[1' 2] 

29 

[1].  A.  Bao,  Ph.D  thesis  (2008)  OSU  [2].  M.  Uddi  et  al,  PCI  32(2009)  929 


Recombination  reactions 

Rate  (cm3s‘1) 

e  +  02(+)  ->  20 

5.6xlO-6T-°5 

He(+)  +  e  +  M->  He  +  M 

1.4xl0-8 

e  +  02  +  M  02(-)  +  M 

4.2xlO-27T-1 

e  +  CH 4(+)  CH3  +  H 

1.0x10-8 

Simulation  results 


A'oa  =  0.34,  XCH4  =  0.16,  P  =  72  Torr,  24  kHz,  a  =  400  1/s 


no  flame,  but  reaction  zone  was  built  up  by  radicals  generated  from  plasma 
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Path  flux  analysis 
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Change  of  branching  ratio 


Change  of  the  branching  ratio  at  the  reaction  zone! 


Reactions 

Normalized 
branching  ratio 

h  +  o2  =  o  +  oh 

1 

e  +  02  =  0  +  0(D)  +  e 

0.48 

e  +  02  =  0  +  0(+)  +  e 

0.42 

e  +  CH4  =  CH3  +  H  +  e 

0.22 

-  1.7 

He(+)  +  02  =  0  +  0(+)  +  He 

0.52 

e  +  02  —  20  +  e 

0.06 

H  +  02  +  M  =  H02  +  M 

0.2 

3 

O 

C 

o 


E 

u 

« 

I 

O 


a  S  curve  transition 


Fuel  mole  fraction 


76%  of  O  production  by  e  and  ions  from  plasma 

Radical  generation  initiated  the  reaction  zone  and  controlled  the  transition!! 


Thrust  3.  Plasma  flame  chemistry  study  in  a  flow  reactor  with 
Molecular  Beam  sampling  Mass  Spectrum  (MBMS) 
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Characteristic  of  low  T  chemistry 
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Transition  from  low  T 
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H202:  low  T  chemistry  indicator 


Difficulties: 

Absorption  -  overlap  with  H20 
GC  -  decomposition/low  reactivity 
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(p=l,  n-heptane/air,  1  atm 
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Molecular  Beam  Mass  Spectrum 


Sampling  system 
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Schematic  of  experiments  with  MBMS 
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Flow  tube  experiments 

DME:  rich  low  temperature  chemistry 
Pressure:  1  atm 


H202  measurement 


Temperature  (K) 


DME  model:  Zhao  et  al.,  Int.  J.  Chem.  Kinet.,  (40)  2008 
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Flow  tube  experiments 
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Conclusions 


1 .  Plasma  can  significantly  accelerate  the  fuel  oxidization  at  low  temperature  to  extend 
the  extinction  limit  dramatically. 

2.  Major  kinetic  pathways  in  plasma  assisted  combustion  were  identified  . 

3.  A  new  counterflow  burner  with  in  situ  discharge  was  developed.  This  burner 
provides  a  new  platform  to  study  kinetic  effect  of  plasma  assisted  combustion. 

4.  The  In  situ  discharge  can  maximize  E/N  at  high  T  flame  region,  therefore,  maximize 
the  electron  energy  and  effect  on  reaction  zone,  and  enhance  ignition  and  extinction. 

5.  The  In  situ  discharge  can  dramatically  enhance  the  ignition  and  modify  the  classical 
S-curve  to  be  a  monotonic  curve. 

6.  MBMS  was  developed  and  H202  was  successfully  measured  directly  for  the  first  time 
in  reacting  system,  enabling  diagnostics  of  intermediate  species  in  plasma  assisted 
combustion  at  low  T. 
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Future  work 


Plasma  part: 

1 .  OH  PLIF  for  counter  flow  diffusion  flame  with  in  situ  discharge  and  compare 
with  simulations 

2.  Low  temperature  plasma  assisted  combustion  for  large  alkanes 

3.  Flow  reactor  experiments  on  liquid  fuel  with  QCL  diagnostics  on  H20,  H202 
and  H02 

4.  Develop  validated  plasma  flame  models 


MBMS  part: 

1 .  Develop  a  JSR  to  study  the  low  temperature  and  high  pressure  chemistry 

2.  Integrate  JSR  with  plasma  discharge  to  investigate  plasma  chemistry 

3.  Develop  advanced  light  source  to  ionize  the  molecular  beam 
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Thanks  the  support  from  AFOSR! 


Questions? 
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Flow  tube  experiments 


H202  calibration 


Dissociation: 
do  it  quickly 

changing  H202  concentrations 
monitor  02  peak 


450  K 


Dilution  gas 
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